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We formulate the low energy limit of QCD like partition functions with bosonic quarks at nonzero 
chemical potential. The partition functions are evaluated in the parameter domain that is dominated 
by the zero momentum modes of the Goldstone fields. We find that partition functions with bosonic 
\Q ' quarks differ structurally from partition functions with fermionic quarks. Contrary to the theory 

with one fermionic flavor, where the partition function in this domain does not depend on the 
chemical potential, a phase transition takes place in the theory with one bosonic flavor when the 
■ chemical potential is equal to m v /2. For a pair of conjugate bosonic flavors the partition function 

shows no phase transition, whereas the fermionic counterpart has a phase transition at [i = m w /2. 
The difference between the bosonic theories and the fermionic ones originates from the convergence 
requirements of bosonic integrals resulting in a noncompact Goldstone manifold and a covariant 
derivative with the commutator replaced by an anti-commutator. For one bosonic flavor the partition 
function is evaluated using a random matrix representation. 



I. INTRODUCTION 



> 

Bosonic quarks appear in QCD whenever the weight of the partition function includes an inverse determinant of 
the Dirac operator. Ratios of determinants occur frequently in analytical approaches to QCD. Inverse determinants 
are for example used to quench a flavor from the theory This form of quenching is useful in describing quenched 

and partially quenched |3|, Li lattice results and is an essential ingredient when computing the spectral properties of 
the QCD Dirac operator. In particular, if we consider the spectral correlation functions of the Dirac operator in the 
""^5 ■ microscopic limit 0, |(| , the theories with bosonic flavors are not only a part of the calculation, they are also part of 
*S ! the result 0,H0. 

The microscopic limit is an extreme version of the chiral limit where the sea and valence quark masses times the 
volume are kept fixed as the volume is taken to infinity. In this limit the zero modes of the Goldstone bosons associated 
. with chiral symmetry breaking dominate the low energy partition function which reduces to a group integral uniquely 
determined by the pattern of chiral symmetry breaking [lOj. Sometimes this limit is referred to as the e-Hmit but 
since the term microscopic limit appeared earlier in the literature (see for example |lllll2j|') we prefer to use this term. 
■ For a review of QCD in the microscopic limit see [l3|. 

In this paper we will consider the low energy limit of QCD at nonzero chemical potential. For fermionic quark flavors 
the baryon chemical potential does not affect the low energy effective partition function at zero temperature. The 
reason is that the lightest degrees of freedom (the pions) have zero baryon charge. Because the fermion determinant 
is complex, QCD at nonzero baryon chemical potential cannot be simulated directly on the lattice by probabilistic 
methods. That is why sometimes a theory is considered where the fermion determinant is replaced by its absolute value. 
For an even number of flavors this corresponds to the product of a fermion determinant and its complex conjugate in 
the weight of the partition function. The low energy partition function for QCD with flavors and conjugate flavors 
depends on the chemical potential p^ |. The conjugate fermionic flavors correspond to ordinary fermionic flavors with 
the opposite sign of the chemical potential |15| . A theory with one fermionic flavor and a conjugate fermionic flavor is 
therefore identical to a theory with two fermionic flavors at nonzero isospin chemical potential. Since the pions have 
nonzero isospin charge this low energy effective partition function depends on the chemical potential. In particular, 
a phase transition to a Bose-Einstein condensed phase takes place at zero temperature when the isospin chemical 
potential reaches m ff /2, see [bil Ha. H71 113. l2Cj . 

In this paper we examine the properties of the low energy QCD partition function with bosonic quarks at nonzero 
chemical potential. We find that the bosonic theories behave completely different from their fermionic counterparts. 
The bosonic partition functions depend on the chemical potential even in the absence of conjugate bosonic flavors. 
By explicit computation of the partition function with one bosonic flavor in the microscopic limit we find that this 
theory has a phase transition when the chemical potential, /i, reaches On the contrary, in the theory with one 

conjugate pair of bosonic quarks we find that the free energy as well as its derivative are continuous at /i = m T /2 
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and zero temperature. The difference in the phase structure between the bosonic and the fermionic theories has its 
origin in the convergence requirements of bosonic integrals which forces us to rewrite the determinants in the partition 
function as the determinant of a Hermitian matrix (also known as hermitization pH I22I EsT]). This procedure leads 
to a noncompact Goldstone manifold and, as will be shown here, a covariant derivative where the commutator is 
replaced by an anti-commutator. A comparison of bosonic and fermionic partition functions is given in table IJ 

In order to compute the partition function with one bosonic flavor we make use of the random matrix representation. 
All other partition functions discussed here have been derived directly from the low energy effective theory. 

The outline of this paper is as follows. We first discuss in general terms the low energy physics of QCD with 
bosonic quarks at nonzero chemical potential. As a warm up exercise we review results for fermionic quarks at 
nonzero chemical potential. We then turn to the theories with conjugate pairs of quarks. In section IVTI we present a 
random matrix model and perform the calculation of the partition function with one bosonic flavor in the microscopic 
limit. Taking the thermodynamic limit of this result then allows us to examine the phase diagram of QCD with one 
bosonic quark. 



II. GENERAL DISCUSSION 



In this section we give a general discussion of QCD-like partition functions at nonzero chemical potential. The 
Euclidean QCD partition function at nonzero chemical potential for Nf flavors with mass m is given by 

Z N f (m; n) = (det^ (D + M 7o + m)) . (1) 

Here and below (. . .) denotes the average with respect to the Yang-Mills action. If Nf > the quarks are fermions 
while for negative Nf they are interpreted as bosons. With fermionic quarks at zero temperature this partition function 
is independent of the chemical potential for [i < mjv/3, which immediately follows from the definition of the grand 
canonical partition function as an average over Boltzmann factors and fugacities. In terms of the representation 
the //-independence of the free energy in the thermodynamic limit can be understood from the gauge transformation 

D + /i 7 o+m = e-^[ J D + m]e' lt , (2) 

so that the factor exp(±/ii) can be absorbed in the boundary conditions of the fermionic fields. In the phase that is 
not sensitive to the boundaries, the partition function does not depend on /i. A nonzero baryon density is obtained 
by baryons winding around the torus in the time direction. In this phase the boundary conditions are important and 
the partition function becomes /^-dependent. 

The second partition function we consider is the phase quenched partition function (the superscript n counts the 
pairs of conjugate fermionic, n > 0, or bosonic, n < 0, determinants) 

Z n=1 (m 7 m*;n) = (det(D + fij + m)det(D + fj,j + m)*) = (dct(D + fij + m)dct(-D + /*7o + m*)), (3) 

so that /1 can be interpreted as the isospin chemical potential Also in this case the chemical potential can 

be gauged into the boundary conditions, this time by a gauge transformation in isospin space. For fi < m n /2 the 
partition function is //-independent at zero temperature. Again, this also follows from the zero temperature limit of 
the Boltzmann factors. 

Next, let us consider the bosonic partition function 

Z (m,m ,fi) - ( det (£ ) + M7o + m ) det ( jD + M7o+m )t' > ~ ^dct(D + /!7o + TO)dct(-Z) + /i7o + TO*)^ ^ 

Because of the nonhermiticity, the inverse determinants cannot be written directly as a convergent bosonic integral. 
However, this can be achieved by introducing the infinitesimal regulator e 

The parameter e may be regarded as source for pions composed of bosonic anti-quarks and conjugate bosonic quarks. 
Expressing the partition function as an integral over the eigenvalues of the Hermitian matrix in (J5J one can easily 
convince oneself that the partition function diverges logarithmically in e. In § this was shown explicitly by performing 
the integrals in the low-energy limit of this partition function. Because of the e-term it is not possible to gauge away 
the //-dependence of the partition function, and we do not expect to find a phase where the partition function is 
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/i- independent. Indeed, we show below that the partition function J5J is /i-dependent for all values of fi. Therefore 
the pion condensate, which has e as source term, is nonzero for all values of /x. In particular this implies that there 
must be a massless mode in the spectrum. We identify this mode explicitly below. 
Finally, we discuss the partition function 

Z N t=~ 1 (m;n) = ( - — i 1 — V (6) 



det(D + fi~fo 

It is not possible to write this partition function as a convergent bosonic quark integral. To properly define the 
partition function we have to rewrite it as 

7 N f =-i( \_/ det(-D + ^70 + m*) 



Z N f = ~ 1 ( m ■ ) — / V P fU ' \ (7) 

( .W- \det(D + /i7o + m)det(-£) + M7Q + m*)/ (> 

so that the inverse determinant can be represented as a bosonic integral after it has been regularized as in JSJ). The 
particle content of this partition function is one conjugate fermionic quark, one bosonic quark, and one conjugate 
bosonic quark. Because of the extra determinant in the numerator, the limit e — > is finite in this case, i.e. the leading 
term is of order e . In the partition function the gauge symmetry therefore allows us to transform the \i dependence 
to the boundaries allowing for a phase transition to occur. As we will show below a phase exists where this partition 
function does not depend on fi. For \x > this phase gives way to a /j dependent phase. 

A summary of the different partition functions discussed in this section is given in the table below. 



Theory 


Number of Charged Goldstone 
Modes for /j, < \i c 


Critical Chemical Potential 


(det(D + /i7o + m)) 





He = \m N 


<|det( J D + /i7 + m)| 2 > 


2 




I 1 ) 

\det(D+^7o+m) / 


4 


He 2^'i' 


( | det(_D+/x7o+m)| 2 ) 


non applicable 


n c = 



TABLE I: Summary of properties of low energy QCD at nonzero chemical potential and zero temperature. 

Notice that the regularization enters because the chemical potential breaks the hermiticity properties of the Dirac 
operator. If we where to consider an imaginary chemical potential, the Dirac spectrum would remain on the imaginary 
axis and no regularization of the bosonic theories would be required [24[ . 



III. FERMIONIC PARTITION FUNCTIONS 



In the microscopic limit the zero momentum modes of the pions dominate the low energy effective partition function 
of QCD in the sense that the nonzero momentum modes factorize from the partition function leaving us with a group 
integral over the zero momentum modes. This integral is uniquely determined by the pattern of chiral symmetry 
breaking, and in the case of bosonic quarks, by the convergence of the integrals. Before turning to the theories with 
bosonic quarks, as a warm-up exercise, we recall results obtained for fermionic quarks. The emergence of the phase 
structure is discussed as well. 
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With fermionic quarks the partition function is automatically convergent and we need only worry about the sym- 
metries. The Lagrangian is determined by local gauge invariance in isospin space ,16]. For Nf + n ordinary fermionic 
quarks and n conjugate fermionic quarks it is given by 

C = ii^TW^C/V^- 1 - ^TrAf (U + U~ l ) (8) 

with U S SU(Nf + 2n), and the quark mass matrix is given by Af=diag(mi, . . . , rriN } , {z} n , { z *}n)- The charge 
matrix B is the diagonal matrix B = diag(l, • • • , ljv/+n: — lj 1 • • > — In)- The chiral Lagrangian is parameterized by 
two low energy constants, the chiral condensate, E, and the pion decay constant F„. The covariant derivatives are 
defined by 

\7„U = d»U + ti5 v0 [U,B}, VvU- 1 =d„U- 1 + t i5 u0 [U- 1 ,B}. (9) 

In the microscopic limit the nonzero momentum modes factorize from the partition function |10| so that the zero 
momentum part of the partition function in the sector of zero topological charge is given by |49| 

Z N f> n ({m f },z,z*;n) = f dU e -^ 2 ^B][u-\B] + ^VTrM(u+u--) _ (10) 

Ju£U(N f +2n) 

Note that the chemical potential and the quark masses only appear through the dimensionless combinations 

fi^iiF^W and M = MT,V, (11) 

where V is the volume of space-time. In the absence of conjugate quarks the charge matrix B is the unit matrix and 
the fermionic partition function is independent of the chemical potential. For example the partition function with 
one fermionic flavor is |{| 

Z N f= l (m) = J (m). (12) 

An explicit expression for fermionic partition functions was derived in |29| . Here we take a closer look at the 
simplest nontrivial example which is given by the theory with one fermionic quark and one conjugate fermionic quark 
both with real mass m, 

Z n=1 (m, m; fi) = 2c 2 ^ / dt t e^* 2 I (mt) 2 . (13) 
The chiral condensate follows by differentiation of the free energy, 

S " =l ^ iA) = IdmlogZ^i^rh^). (14) 

In figure^we have plotted the chiral condensate as a function of m (in units of m c — 2[i 2 F 2 /£) and as a function of 
fi (in units of /i c = \J mE / \\[2F 1T )') . The full curves display a smooth dependence as was to be expected for a finite 
volume. In the thermodynamic limit, the would-be phase transition at m = m c in the left figure or at /i = fi c in the 
right figure becomes sharp. In terms of the dimensionless variables this limit is reached for rh — > oo or /t — > oo and a 
kink develops at the expected value of m/(2/i 2 ) = 1. 

In order to see how the kink is recovered we take the thermodynamic limit of l|13(l which is given by a leading order 
saddle point approximation. Using the asymptotic form of the Bessel function we obtain 

Z n=1 (m,m;/i) -e 2 ^ / dt — e^* 2 e 2t ™ . (15) 
Jo m 

A phase transition occurs when the saddle point, 

< = W m 

hits the integration boundary at t = 1, that is when rh = 2/t 2 . In the phase where the saddle point is inside the 
boundary we obtain 



]™ _1 (m;/i) m 
E 24l 



(17) 
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If the saddle point is outside the integration domain the chiral condensate is simply given by 



(18) 



These results are shown by the dashed curves in figure Q and are in agreement with the results obtained from chiral 
perturbation theory [Tl fl7l 111 IT^ . 




FIG. 1: The chiral condensate for n = 1 (one fermion and one conjugate fermion). Left: The mass dependence for fi = v5 
is shown by the full curve. In the thermodynamic limit (ft, — > oo, dashed curve) the condensate grows linearly with m until it 
reaches a plateau at m = m c = 2/i 2 i^/E. Right: The full curve displays the chemical potential dependence for m = 10 and 
the dashed curve shows the result in the thermodynamic limit (m — ► oo). 



IV. THE THEORY WITH A PAIR OF CONJUGATE BOSONIC QUARKS 

We now turn to QCD with bosonic quarks at nonzero chemical potential and start off with the theory with a pair 
of conjugate bosonic flavors JIJ. This theory is simpler than the theory with a single bosonic quark since the latter 
involves both fermionic and bosonic flavors (see 10). 



A. Covariant Derivatives 



As discussed in section [H] the n = — 1 partition function diverges. To derive the effective partition function we 
therefore need to take into account the convergence of the integrals in addition to the symmetries. If we write the 
Dirac operator as 



D = 



d 
-df 



(19) 



the determinant in © can be rearranged as 



e 
e 

z* —d + /i 
d< + /! z* 



d + ix 

(1 z 



e 



e z d + [i 

z* e d — (i 

(v + /i e — z* 

$ — [a —2 e 



(20) 
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For the purpose of studying the transformation properties, we rewrite the determinant of the right hand side in a 
2x2 block notation as 



Ci d + Bi 
d) + B 2 C2 



(21) 



This operator becomes locally gauge invariant under time dependent but spatially constant flavor gauge transforma- 
tions 



Ci d + B l 
d^ + B 2 C2 




if B\ and B2 are transformed as 



and the mass matrices are transformed as 



(22) 



Bx ->■ vBxv^ 1 - [dovjv^ 1 

B 2 -> UB2U- 1 + \dou\u- 1 (23) 



Ci -> vC,iu 1 

C2-»uC2« _1 . (24) 

In Q we showed that the Goldstone manifold is this case is given by Gl(2)/U(2). The transformation (|22|l induces 
the following transformation on the Goldstone fields 

Q^uQv- 1 . (25) 

Therefore, 

doQ — > OquQv^ 1 + udoQv -1 + uQdov^ 1 , 
doQ- 1 -> dowQ" 1 ^ 1 + waoO" 1 ^" 1 + «Q _1 a M _1 . (26) 

One immediately sees that the covariant combinations are 

d Q - QBx - B 2 Q, 
d Q- 1 +B 1 Q- 1 + Q- 1 B 2 . (27) 

Of course B\ v — B 2 ^ — [ib V Qaz = fxS^oB. We thus obtain the covariant derivatives 

V v Q = d v Q-fi5 v o{Q,B}, X7.Q- 1 =d u Q- 1 +fi5 v0 {Q- 1 ,B}. (28) 

The chiral Lagrangian is the low-energy limit of QCD and should have the same covariance properties as were 
derived above. Taking into account terms to order p 2 in momentum counting we find the Lagrangian 

F 2 i 
C = -— TrV^QV.Q- 1 - -STrM (Q - IQ- 1 !) (29) 



with 



z e 




M = M Z ) and /=("*). (30) 



As before, the pion decay constant is denoted by and the absolute value of the chiral condensate is given by S. We 
emphasize that the sign before the kinetic term is opposite to the sign of the kinetic term in the fcrmionic Lagrangian. 
This sign enters to compensate the sign due to the noncompactness of the modes in Gl(2)/U(2) so that the kinetic 
terms pion fields have the correct sign [2{| [2(J . 

In the case of the fermionic partition function with one flavor and one conjugate flavor, we also could have rearranged 
the fermion determinant as in (|20[1 . and we would have obtained covariant derivatives as in l)28[l. Indeed, if we make 
the transformation 

U^Uax. (31) 
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in the Lagrangian JSJl, the covariant derivatives in © change to 



d v U + iiS v0 [U,B] {d„U - ^vo{U,B})ax, 
d v XJ- x +{jl8 u0 [U-\B] a^O.U- 1 +n8 v0 {U- 1 ,B}). (32) 

Notice that both U and Uo\ are unitary. The mass term in (jHJ becomes 

Tr ^ °^ (C/ + U- 1 ) = TrM(U + a 3 U-'a 3 ) - Tr f * £ ^ (C7 - ^E/" 1 ^)^ = Tr (J Z J ^ ([/ + It/- 1 /). 

(33) 

The invariance properties of this term are consistent with those of the representation l|2U|) . Because both U and — U 
belong to U(2) this term has the discrete symmetry M — > —I MI. In the bosonic case this symmetry dictates that 
the relative sign of the two terms in the mass term has to be negative. 



B. Microscopic limit of the partition function 

In the microscopic limit, the partition function is dominated by the zero momentum modes and in the sector of 
zero topological charge it is given by 

Z n =-\z,z*;n) = lim / -^-9(Q) e -%Kn 2 Tr{Q,B}{Q-\B} + fvinw^CQ-rtr^ (34) 
e->oJ det Q 

where dQ6(Q) /det 2 Q is the integration measure on positive definite 2x2 Hermitian matrices [5(|. In the limit e — > 
the integrals can be performed analytically resulting in [j| 

Z—^i*;/*) =loge^ex P Ko (jf ) . (35) 

Using the Toda lattice equation, the product of this partition function and its fermionic counterpart gives the quenched 
microscopic spectral density at nonzero chemical potential. This has been confirmed beautifully by recent quenched 
lattice QCD simulations at nonzero chemical potential with a staggered Dirac operator |2?| and an overlap Dirac 
operator psj . 

The volume dependence has dropped from the partition function l|35(l and no kink will develop in the thermodynamic 
limit. The chiral condensate is given by 

v — x (m;/i) _ 1, , . /: , : 



-drh log Z n (rn, m; /}) 

rh rh 
~ ■— r for ui — > oo , fi — > oo and — ^ ~ 1. (36) 
2/x 2 2/i z 

It follows that the low energy theory with one bosonic flavor and its conjugate does not have a phase transition as 
function of ft. This is to be contrasted with the fermionic counterpart which, as discussed in section ITTT1 has a phase 
transition at rh = 2/2 2 . Plots of the chiral condensate for n = — 1 are shown in Fig. 

The reason for the absence of a phase transition at rh = 2/t 2 becomes clear from the calculation of the mass spectrum 
of the Lagrangian (|29J) along the lines of ^3] (see Appendix A). For all values of the chemical potential we find a 
charged massless mode (in fact with a mass ~ \fe) which condenses for any [i > 0. The masses of the remaining three 
modes are given by lix . rh j/2ii and 2/x-\/l + 3(m 7r /2/i) 4 , exactly the same as in other QCD-like theories at nonzero 
chemical potential |17l llflj . 



V. THE EFFECTIVE PARTITION FUNCTION AND CHIRAL RANDOM MATRIX THEORY 

We now turn to the computation of the partition function with one bosonic quark, that is Nf = — 1. As follows from 
eq. Q the Goldstone bosons result from the breaking of chiral symmetry in the theory with one fermionic and two 
bosonic quarks. The partition function is therefore given by an integral over the super group G/(l|2). Such integrals 
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FIG. 2: The chiral condensate for n = — 1 (one boson and one conjugate boson). Left: The mass dependence for p, = \fE 
(full curve). Right: The chemical potential dependence for rh = 10 (full curve). In this case there is no phase transition 
at (i = nij/2 as can be seen from the result in the thermodynamic limit (dashed curves, obtained for /} — *■ oo or m — > oo, 
respectively) . 



are technically rather complicated, and at this point we have not succeeded to evaluate the Nf — — 1 partition function 
along these lines. 

Fortunately, the effective low energy partition function has an alternative representation as the large N limit of a 
random matrix theory with the same chiral symmetries. This has been shown explicitly in the fermionic case at zero 
and nonzero 0, 123 . 1.SCj chemical potential as well as for bosonic |3lj| and supersymmetric [82l l3.il l35| partition 
functions at zero chemical potential. For bosonic quarks at nonzero chemical potential the equivalence between the 
effective theories and random matrix theory in the microscopic limit has been established for n = — 1 in |36| . 

As before, we only consider the theory in the sector of zero topological charge. 



VI. THE RANDOM MATRIX MODEL 



The random matrix partition function with Nf quark flavors of mass m and n pairs of regular and conjugate quarks 
with masses y and z* , respectively, is defined by 30] 

ZN ,n ({ m f}>V> z *->ti = /«^(^GWdet^(%)+ m/ ) 

x det" (V(fi)+y )det n ( Z>t ( M ) + **) , (37) 
where the non-Hermitian Dirac operator is given by 

Here, $ and iff are complex N x N matrices with the same Gaussian weight function 

w G (X) = exp(-A^TrX t X) . (39) 

Bosonic quarks appear as inverse determinants and notationally we simply allow Nf and n to take negative values. 

Of course, in general the QCD partition function and the random matrix partition function are different. However, 
when we consider the microscopic limit where the variables 

m = 2mN and p 2 = 2fj, 2 N (40) 



8 



are fixed as N — > oo the random matrix partition function and the QCD partition function coincide provided that we 
identify (see the discussion in (3|J) 

TO = 2mN -> mVY> (41) 
fi = 2^i 2 N -> fi 2 F 2 V. 

In this section we will work within the random matrix framework and use the identifications (|41() in the final results. 
In |30l ] it was shown that the random matrix partition function (|37|) can be rewritten in the eigenvalue representation, 

r N 

Z% f ' n (m,y,z*;n)~ / J[ d 2 z k V N ^ n ({z t }, {z*}; M ), (42) 
Jc k=i 

where the integration extends over the full complex plane and the joint probability distribution of the eigenvalues is 
given by 

1 N 

V N f> n {{ Zi },{zt}^) = —lA^azf})! 2 l[w(z k ,z* k ;n)(m 2 - z 2 k ) N t(y 2 - zl) n (z* 2 - z* k 2 ) n . (43) 
^ fe=i 

The Vandermonde determinant is defined as 

N 

A N ({zf}) = J] (z 2 - z 2 ), (44) 
i>j=l 

and the weight function reads 

w(z k ,z k ;(i) = \z k \ K I ^ |Zfc| I cxp I — (z fc + z fc )l. (45) 

The eig envalue re p resentation makes it possible to employ the method of orthogonal polynomials in the complex plane 
|3^. i38L I39L I40L l4ll| to compute the spectral density and eigenvalue correlation functions |3Cj . 



A. Orthogonal Polynomials and their Cauchy transform 

In order to evaluate the partition function with Nf = — 1 we will make use of orthogonal polynomials and their 
Cauchy transform. The complex Laguerre polynomials given by [3(| 

p k (z;n) = (^-) k\L k (-r^s) (46) 



N J \ 1- n 

are the orthogonal polynomials corresponding to the weight w(z, z*\ fx) given in (|45|l . To be specific, the polynomials 
satisfy the orthogonality condition 42] 

<Pz w(z,z*;n) p k (z;n) Pi(z; /j,)* = r k 5 k i , (47) 



with the norm 



ir^ (l + ^ k k\ fc! 
rk = JpkT2 • ( 48 ) 

The Cauchy transform of the orthogonal polynomials is defined as 

hk(m;n)= / d 2 z-^— — ~w{z, z*; n)p%{z; fi), (49) 
J c z z — 

where C indicates that the integration extends over the complex plane. Using that the weight function and polynomials 
are even functions of z the Cauchy transform can be written as 

h k (m;n)= [ d 2 z——^- riv(z, z*; fi)p* k (z; (1). (50) 

Jc z(z-m) 
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VII. QCD WITH ONE BOSONIC FLAVOR 



It was shown in 361 that 



Z 



N f = -1 
N 



(m;/x) 



TN-l 



W-l 



(m;/i). 



(51) 



Therefore, studying the properties of bosonic partition functions is equivalent to analyzing the properties of the 
Cauchy transform. While the relation between the partition function and the Cauchy transform was established in 
[36| . no explicit evaluation of the Cauchy transform was given; this evaluation follows below. We will find that the 
partition function with one bosonic quark depends on the chemical potential for /i > m v /2. 

We are interested in the microscopic limit where N ~ > oo for fixed rh — 2Nm and fixed /t 2 = 2[i 2 N . (This scaling 
of \x is the analogue of the weak non-hermiticity limit introduced in |43j.) In this limit the polynomial pjv— i/riv— l is 
given by 



pjy_i(z;/i) 



1 (N-l)\ fi\ N -i r 

~ {1 ~2N> Ln - 1 



rjv-i N N 
2N 5 l 2 e N 



AN J 



where we have used that 



/2irirfl 2 ' 
ttA 2 (1 + f^Nf^-^iN - 1)\{N - 1)! 



2N 2N+1 

The microscopic limit of the partition function is thus given by 



Z N ?=~ 1 (m;fL) = - 



2N l / 2 e N 

^2 ' 



d 2 z- 



To calculate this integral we write (with z — x + iy) 



I (z") - Sgn(y)-(K (z*) - K (-z*)) 



(52) 



(53) 



(54) 



(55) 



and first calculate the integral over y by a complex contour integration. The first term in (|55|l is exponentially damped 
in the upper half of the complex y-plane and the second term in the lower half-plane. This allows us to close the 
integration contour of the y variable, and by Cauchy's theorem we obtain for the real part of the partition function 



Re(Z N f— 1 (m;iJ.)) = 2nc N (fi) dx(2x-m)exp 



2x 2 — 2xrh + m 2 



4fi 



;,2 



(rh — x)6(2x — m)Ko 
-6{m - x)6{m - 2x)I 
The imaginary part of the partition function is given by 
l-m{Z N i=- 1 {m-jl)) = 

(x — m)K 



2cjv(/t) / dx(2x — rh) exp 
(2x — rh)rh 



(2x — rh)rh 
/ (2x — m)rh 

{ w 

2x 2 — 2xrh + rh 2 



Io{2x — rh) 
Ko(rh — 2x) 



(56) 



Afi 2 

+9(rh - x)B(2x - rh)K 
+9(rh - x)6{fh - 2x)K 



4(l 2 

Kq(2x — rh) 



(2x — rh)rh 

(rh — 2x)rh 
4^2 



Kq(2x — rh) 
KQ(rh — 2x) 



(57) 
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The overall constant cjv(/t) is defined by 



cat (A) 



2N 1 / 2 e N 

r 2 ' 



(58) 



The integrand of the imaginary part is odd about x = to/2 so that the integral vanishes. Using s = x — to/2 as 
new integration variable we obtain the final expression for the partition function with Nf = — 1 at nonzero chemical 
potential 



m/2 



dssexp[-^} 



(«)*> ( ^ ) Io(2s) 



MA,(|p)/v 1 ,(-2.) 



47TCAr(/i)e s f' 2 



m/2 



c 

dssexp[-^^]X 



STO 



/o(2s) 



dssexp[-— ]/ 



sm 



) K (2s) 



(59) 



A. The chiral condensate for one bosonic flavor 



The chiral condensate given by 

£Af/ ^ ( " ;A) = -a m log z»'-H*i A) (60) 

is plotted in figure (full curve). Notice that the chiral condensate depends on the chemical potential, and that for 
m 3> 2/t 2 , it approaches the value S. Below we will show that in the thermodynamic limit E^^" 1 develops a kink 
at to = 2/i 2 . 

1-5 I 1 1 1 1 1 1 1 1 1.25 , . , . , . , . , . , 




FIG. 3: The chiral condensate in the theory with one bosonic flavor. Left: the mass dependence for ft — v5 (full curve). In 
the thermodynamic limit (p, — > oo, dashed curve) E -1 changes from a linear dependence on m to become m independent 
for m > m c = 2/x 2 i^/E. Right: The dependence on the chemical potential for rh = 10 (full curve) and in the thermodynamic 
limit (rh — > oo, dashed curve). 



B. Phase transition with one bosonic flavor 



Above we have determined the partition function in the microscopic limit, where the source m times the volume 
N, is kept fixed, and as a consequence, the phase transition at rh = 2/i 2 is smeared. A cusp in the derivative of the 
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free energy only appears in the thermodynamic limit where N — > oo for fixed m. This limit can also be approached 
by taking m = 2mN —>■ oo and fi 2 = 2/i 2 A — > oo in the microscopic results as will be done below. 

In the thermodynamic limit the integral (|59|l can be calculated by a saddle point approximation. The saddle point 
of the first integral is at 

s = -- + 2fi 2 , (61) 

and the saddle point of the second integral is at 

-2 ™ 



-2/i 2 + -. (62) 



In the thermodynamic limit the partition function for rh > 2/i is therefore given by 

Z N f=- 1 {rh;fi) = ^LcjvOjVV^™. (63) 

The exponent depends linearly on rh resulting in a chiral condensate that is equal to 1. In the normalization for 
which the partition function is /i-independent for small \i (see pMj), the exponent in c^iji) (see eq. (I58|H becomes 
exp(— 2/i 2 ) instead of exp(— |/t 2 ). A phase transition occurs when the dominant saddle point hits the boundary of 
the integration domain, i.e. at rh = 2/t 2 , cf. 1)61(1. For rh < 2/i 2 the saddle point is outside of the domain of the first 
integral and the main contribution to the first integral comes from the region close to rh/2 so that 

Z^=- 1 (m;A)~e" 2 / 2 ' i2 - 2 ^ for fi 2 > y. (64) 
The condensate in the thermodynamic limit is therefore given by 



Nf=-i(^,. r,\ \ _1 _i_ ™_ r^r ™ ^ or,2 
1 for rh > 2fi 2 



E JV /=- 1 (m;/t) I -1 + |l for rh, < 2/i 



(65) 



In figure |3J this result is shown by the dashed curves. 

As promised, we have shown that in the thermodynamic limit the chiral condensate has a kink at rh — 2fi 2 . Using 
(|41ll and the GOR-relation (2mS = F 2 m 2 ) we see that the phase transition takes place at /i = m n /2 with 

S^^'Kiri = f -1 + 5^ for < 2[i ^ 
S 1 for m v > 2/i 

This behavior of the chiral condensate can be explained from the structure in (J7J). We first note that the average mass 
derivatives of fermionic and bosonic determinants give factors that are equal in magnitude but have opposite sign. 
At /i = the contribution from the fermionic quark thus cancels against the contribution from one of the bosonic 
quarks resulting in a total chiral condensate of magnitude E. For /i < m v /2 the partition function is /i independent 
and the chiral condensate therefore remains at the value E. As fj, exceeds m^/2 it becomes energetically favorable to 
create pions made out a bosonic anti-quark and a conjugate bosonic quark. As a consequence the contribution to the 
chiral condensate from the bosonic quarks starts rotating to zero for fi > m^/2 so that the total chiral condensate for 
large values of ix/m^ is given by the contribution of the fermionic flavor alone. For /z ^> m^/2 the chiral condensate 
therefore is equal to the chiral condensate at /i = but with the opposite sign. 
The second order phase transition can also be seen in the baryon density given by 

= -a A logZ (m; M ) = | m £ } for A2> | • (67) 

In contrast, the free energy for Nf — 1 does not depend on jl. 

When rh ~ 1 the finite volume effects for the chiral condensate are strong and the microscopic prediction differs 
significantly from the result in the thermodynamic limit. For m<g 1 the partition function is given by 

Z N s=~ x (rh; fi) = A+ ^ log rh + 0(rh 2 ) (68) 
8 

resulting in a chiral condensate 

E Ar /=- 1 (m;/i) 1 „ , . , 

— = — -mlogm lor m«l, (69 

£ AA v ' 

where A is a constant. This explains the approach to zero for rh — > of the solid curve in the left panel of Fig. [21 
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VIII. THE PHASE STRUCTURE AND THE EIGENVALUE SPECTRUM 



The chiral phases of QCD reflect themselves in the eigenvalue spectrum of the Dirac operator. For zero chemical 
potential the eigenvalue density near the origin of the imaginary axis is an order parameter for spontaneous breaking 
of chiral symmetry |44j | . When the chemical potential is nonzero the eigenvalues move into the complex plane and 
the chiral condensate in full QCD is then linked to the oscillations of the eigenvalue density . For nonzero isospin 
chemical potential (n — 1) the transition into the Bose-Einstein condensed phase takes place when the quark mass is 
at the boundary of the support of the eigenvalue density (see for example |lSj). 

Also in the case with one pair of conjugate bosons, n = — 1, the phase structure is closely related to the Dirac 
spectrum. Setting n = — 1 in (|42|l we obtain the following eigenvalue representation of the partition function, 

1 f N 2 N 1 
z N f -0,n— i( 2)Z * ;M )^ Y[d 2 z k \A N ({zf})\ 2 Y[w(z k ,z* k ;fi)— - 2 -2-. (70) 

A* J c k=1 k=l \ z z k)\ z z k ) 

We observe that the probability to find an eigenvalue, Zj, has a modulus squared pole at the mass z. This gives rise 
to the logarithmic divergence discussed in section II, 

L»/ u v^~w loge - (71) 

Here, C e (u) is an annulus centered at u with outer radius 1 and inner radius e < 1. Because of the Vandermonde 
determinant, the probability that two or more eigenvalues are close to the mass z does not diverge. Therefore the 
logarithmically divergent contribution to the partition function is from eigenvalue configurations with exactly one 
eigenvalue close to the mass. Setting zjq = z in l|70|) and using that 

lAa^}^)! 2 !^^^^ 2 }^ 1 )! 2 !! 1 !^-- 2 ! 2 , (72) 

i=i 

we find that the product on the right hand side exactly cancels the two bosonic determinants. The remaining integral 
over N — 1 eigenvalues is equal to the quenched partition function which does not depend on the mass. The only 
mass dependence is from the weight function, w, evaluated at the mass z and from the divergence (|71|l . We thus find 

which is exactly the result l|35fl given in section IV. 

The absolute value squared pole responsible for attracting the eigenvalue zn to the mass z appears for all nonzero 
values of /i. Hence the n = — 1 theory must be regulated for all nonzero values of [i. The transition to fx — > is 
discontinuous because if jjL = the probability to find an eigenvalue near z is zero unless z is on the imaginary axis. 



IX. SUMMARY AND DISCUSSION 



We have computed the microscopic limit (also known as the e-limit) of the QCD partition function at nonzero 
chemical potential with one bosonic quark. For the computation we made use of the random matrix representation 
of the partition function. Contrary to the partition function with one fcrmionic flavor this partition function depends 
on the chemical potential. In particular it has been shown that a phase transition takes place at ji = The 
theory with a pair of conjugate bosonic flavors also behaves differently from its fermionic counterpart. While a phase 
transition signaled by Bose-Einstein condensation of pions takes place in the fermionic theory, its bosonic variant 
always remains in a \x dependent phase. The main differences between fermionic and bosonic partition functions have 
been summarized in table 1. 

Bosonic partition functions as studied in this paper also appear in the closely related Hatano- Nelson model |4fij . 
This is a model for a disordered system in an imaginary vector potential. The random matrix limit of this model 
is a non-Hermitian random matrix model just like that in l|37H except for the chiral block structure in (|38fl . The 
partition function with Nf = — 1 can again be expressed as the Cauchy transform of the orthogonal polynomials. 
However, because this model is technically simpler than the random matrix model for the QCD partition function, 
this partition function can also be evaluated directly by means of standard random matrix techniques. The agreement 
of both approaches |47) confirms our understanding of the Cauchy transform representation of the partition function. 
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Moreover, as in the chiral case, we find that the partition function of the Hatano-Nelson model with Nt — — 1 depends 
on the chemical potential and a phase transition appears in the thermodynamic limit. 

These results generalize to theories with more fermionic and bosonic flavors, and we expect similar surprises for 2 
color QCD and for QCD with quarks in the adjoint representation. 

Our results might be relevant for a better understanding of the sign problem in QCD at nonzero baryon chemical 
potential. The expectation value of the phase of the fermion determinant is given by 



It thus corresponds to a partition function with a bosonic quark and its conjugate and two fermionic quarks. If the 
same arguments apply as in the case of the partition function with one bosonic flavor the vacuum expectation value 
of exp(2i6>) will be nonanalytic at [i = m n /2, in the quenched as well as in the unquenched theory. We speculate 
that this nonanalyticity persist at nonzero temperatures and that beyond this point the severity of the sign problem 
increases drastically. 

Note added in proof: The expectation that (exp(2i#)) is nonanalytic at /i = m n /2, in the quenched as well as in 
the unquenched theory has now been verified 48]. 
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Appendix A. CALCULATION OF THE PARTICLE SPECTRUM OF THE CHIRAL LAGRANGIAN 

FOR A PAIR OF BOSONIC QUARKS 

In this appendix we calculate the particle spectrum of the chiral Lagrangian Ij29(l for QCD with a pair of conjugate 
bosonic quarks. We only work out in detail the case that z = x is real and e = 0. 
The quark fields in the Goldstone manifold Gl(2)/U(2) can be parameterized as 

q_ l coshse r+t sinhse l9+ * j — ( coa hse~ r ~ t — sinhse* 9- ' \ 

\ smhse _l6,+ * coshse~ r+t I ' I — sinhse~ j6, ~* coshse r ~ t J 

For the covariant derivatives we find 

V Q = d Q-fi{Q,B} 

_ I e r+t [d cosh s + cosh sdo(r + t) — 2/j, cosh s] e i9+t [do sinhs + smhsdo(id + t)} 



\ e l0+t [do sinhs + sinhs<9o(— i0 + t)] e r+t [do cosh s + cosh sdo(—r + t) + 2/i cosh s] 

and 

VqQ- 1 = doQ- 1 +fi{Q-\B} 

_ I e~~ r ~ t {do cosh s — cosh sdo(r + t) + 2/i cosh s] — e ie ~*[9o sinhs + sinhs9o(*# — t)] 



(76) 



(77) 



\ — e 10 t [do sinhs — s'mhsdo(i9 + t)] e r 1 [do cosh s + cosh s8q (r — t) — 2[i cosh s] 
For the time derivatives in the kinetic term we obtain 

TrVoQVoQ -1 = -2(<9 s) 2 - 2(d Q t) 2 - 2 cosh 2 s(d Q r) 2 - 2 sinh 2 s(<9 <9) 2 + 8/i cosh 2 sd Q r - 8fi 2 cosh 2 s. (78) 
The mass term is given by 

— i—T,TrM T (Q — JQ -1 /) = — iS- (Ax sinhs cosh t cos 9 — Ay sinhs sinhi sin 9 + 4e cosh r cosh s cosh t) (79) 
The total Lagrangian is 

F 2 1 

-— TrV^gVQt _ j-£TrM T (Q - IQ" 1 /) (80) 

with M and J defined in (|3(J|) . Notice that because of the noncompact degrees of freedom, the sign of the kinetic term 
is opposite to the usual sign. 

The static part of the Lagrangian for y = and e = is given by 

^static = 2F 2 /i 2 cosh 2 s — i2Hx sinhs cosh t cos 9 — 4e cosh r cosh s cosh t. (81) 

The minimum is at (note that we can always divide out coshs, so there is only one minimum) 

sinhs = =i^| (82) 

and 

< = 0, r = 0, (9 = 0. (83) 

We expand the chiral Lagrangian H80JI to second order in the small fluctuations about the saddle point. The second 
order terms are given by 

F2 [-^(d Ss) 2 - ^(d St) 2 - i cosh 2 s(d 5r) 2 - i sinh 2 s(d S9) 2 + ifi sinhs cosh sSsd 5r - 2 / u 2 (cosh 2 s + sinh 2 s)(Ss) 2 } 
+iZxsmh s[(5s) 2 + (St) 2 - (S9) 2 }. (84) 
We obtain the uncoupled masses 

2 2xE m% 

m t = -z_smhs = ^, (85) 

m 2 = 4/A (86) 
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The remaining two masses of follow from solving 

det 



E — 4/i cosh s +4/ii? sinh s cosh s 
—4(j,E sinh s cosh s E 2 cosh 2 s 



since cosh s ^ we can divide by cosh 2 s and obtain 

E 2 (E 2 - 4^ 2 cosh 2 s) + 16fi 2 E 2 sinh 2 s = (87) 

with solutions 

£ 2 = or £ 2 = 4/i 2 cosh 2 s - 16/i 2 sinh 2 s = 4/i 2 (1 + ^) (88) 

16/i 4 

Note that one mode is massless for all fi > 0. For e ^ one can show along the same lines that the massless mode 
obtains a mass ~ \fi. 
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Without loss of generality we consider the partition function in the trivial topological sector. 

In |3?j| we incorrectly used commutators in I34H instead of anti-commutators. This resulted in an extra overall factor 
of exp{2VF^ 2 ) which was corrected for by hand in the calculation of the spectral density. In Q, the correct form of 
Z n= ~ 1 (z, z*; fi) was obtained because the overall p-dependent constant was obtained by matching to the spectral density 
via the Toda lattice equation. 



